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Abstract. We argue that an ionic lattice surrounded by a Fermi liquid changes phase several times un- 
der pressure, oscillating between the symmetric phase and a low-symmetry dimerized structure, as a 
consequence of Friedel oscillations in the pair potential. Phase oscillations explain the tendency towards 
dimerization which has been recently reported for the light alkali metals under high pressure. Moreover, a 
restoring of the symmetric phase is predicted for such elements at an even higher density. 

PACS. 71.20.Dg Alkali and alkaline earth metals - 62.50.+p High-pressure and shock-wave effects in 
solids and liquids - 64. 70. Kb Solid-solid transitions 



1 Introduction 

Some years ago we predicted £Q that under high pres- 
sure the light alkali metals should be unstable towards 
less symmetric phases. We also argued that the lowering 
of symmetry could give rise to a metal-insulator transi- 
tion. We did not attempt to predict the structure of the 
high pressure phase, and unfortunately any such theoreti- 
cal prediction is quite difficult, because of the small energy 
differences which separate most phases. Our result fol- 
lowed the early suggestion that electron correlation in the 
alkalis could give rise to broken symmetry of the charge 
density wave (CDW) type, albeit incommensurate with 
the lattice [2]. The role of core repulsion in stabilizing the 
high pressure phases of the alkali had been also empha- 
sized by McMahan et al. [HIFT]. 

The prediction of a symmetry lowering was quite un- 
expected, as it is usually believed that an increase of pres- 
sure should give rise to an increase in the symmetry of 
the system. Actually, under very high pressure the most 
likely state of matter is a uniform plasma, so that any sys- 
tem is expected to raise its symmetry at a sufficiently high 
pressure. However, many "simple" metals are still far from 
such a limiting behaviour, and their path towards the ulti- 
mate metallization may contain oscillations of their degree 
of symmetry. 

More recently, both theoretical calculations [3] and ex- 
perimental findings |B]|7|E] have confirmed the tendency of 
the light alkali metals to lower their symmetry under high 
pressure. Low symmetry phases have been observed for 
Li |7] and Na (see also Ref. ^HDi whereas the even- 
tual occurrence of a metal-insulator transition has not 
been observed yet [Sj- According to recent calculations 
both Li and Na show a tendency towards the forma- 



tion of atomic pairs, and a dimerized oC8 structure would 
be the most stable phase above 165 GPa for Li, and above 
220 GPa for Na. 

The physical reason for dimerization is not evident. 
Ab initio electronic structure calculations 12 indicated 
a tendency towards "distance alternation" , due to a size- 
able overlap of pn orbitals in the interstitial regions. On 
the other hand, it has been noticed ^l] that the increase 
of s-p hybridization could give rise to a low coordina- 
tion number. However, even the fully dimerized phase is 
far from any standard covalent solid: the electron density 
is uniformly spread and almost constant, while the first 
and second neighbours distances are comparable. Such a 
dimerized phase is better described as a charge density 
wave in a high density metal rather than a molecular solid. 
Now, at zero pressure both Li and Na are already well 
described by a simple degenerate Fermi liquid where the 
kinetic energy is the most relevant energy term |13| . At 
high density all the interactions should become smaller 
and smaller compared to the kinetic energy, and if all the 
energy terms are assumed to be monotonic functions of 
density, then nothing relevant is expected to happen. On 
the other hand, the high density limit is where the Fermi 
liquid model should work better, thus it remains to be ex- 
plained why a Fermi liquid should be unstable towards a 
charge density wave. 

In this paper we address such a problem and show that 
an ionic lattice, surrounded by a Fermi liquid, oscillates 
between a simple metal and a charge density wave, sev- 
eral times with increasing density. We start in Sec. [21 by 
reviewing the evaluation of the structural energy of a solid 
lattice within linear response theory. In Sec. we later 
discuss the case of lithium and other alkali metals under 
high pressure, where the presence of Friedel oscillations in 
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the screened pair interactions can justify the occurrence 
of an instability towards a broken-symmetry phase. We 
eventually summarize in Sec. fj] 



2 Structural energy of solids close to a phase 
instability 

The basic idea is foreshadowed in the seminal works of 
Pettifor et al. |14U15j , who discussed the possible existence 
of phase oscillations. The physical motivation is related to 
Friedel oscillations rjHj, which characterize the screened 
interactions and give rise to a non-monotonic oscillating 
behaviour of the structural energy terms. A dimerization 
of the lattice can be regarded as a broken-symmetry phase 
described by an order parameter y which vanishes in the 
symmetric phase. In most cases, y can be taken to be the 
difference between first and second neighbours distances, 
and y = yields the monoatomic lattice. The total struc- 
tural energy of the system must be an even function of 
y and can be expanded as U(y) « U(0) + ^U"(0)y 2 . Its 
main contribution comes from screened pair interactions 
which oscillate with a wavevector q = 2kp, where kp is 
the Fermi wavevector (Friedel oscillations). Thus the sign 
of the second derivative U" (0) is dominated by the sign of 
the second derivative of the pair interaction with respect 
to the pair distance. But an oscillating pair interaction 
yields oscillating derivatives, so that U"(0) is expected to 
change sign as the pair distance decreases. A positive sign 
U"(0) > corresponds to having a minimum for y = 
(symmetric phase); a negative sign U"(0) < corresponds 
to having a relative maximum for y — 0, thus indicating 
that the symmetric phase is unstable towards a dimerizcd 
phase. Moreover, we expect that with increasing density 
the sign of U"(0) could change several times, thus giv- 
ing rise to phase oscillations between a symmetric and a 
dimerized lattice. 

In order to make the above idea more quantitative, 
let us specialize to the light alkali metals. We start by 
expressing the total energy per atom of a metal lattice in 
real-space formulation as |17j : 



(i) 



where the pair potential <P(Rij) measures the interaction 
between two ions located at sites i and j in the lattice, i?^ 
being their mutual distance, and Uq summarizes all the 
contributions independent of the lattice structure. Within 
second order local pseudopotential theory, the pair poten- 
tial can be written as |15U18| : 
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R 



dq 



where Z is the atomic number, 
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with e(q) the dielectric screening function, n 2 — Akp/n is 
the Thomas-Fermi screening parameter, x{o) is the nor- 
malized Lindhard susceptibility [x(0) = 1], and G(q) takes 
into account for exchange and correlation corrections (lo- 
cal field corrections) to the electron-electron interaction 
|19j . In the following, we shall take Z = 1 for simplicity, 
which is a good approximation for the light alkali at am- 
bient conditions. In Eq. J2J) we assume the 'empty core' 
model for the ionic pseudopotential |20| . 



v(q) = cos qR c 



(4) 



where R c defines the radius of the atomic core. The core 
radius R c is usually obtained by fitting Eq. fJl against 
the value of the band gap 2|u(g)|. We have checked that 
Eq. (TJJ is also consistent with the calculated band struc- 
ture of lithium both in the bec |21| and in the fee struc- 
ture |22| . However, different sources of experimental data 
result in slightly different values of R c for a given element 
|23j . Moreover, the value of R c can be effectively altered 
by chemical substitution (e.g. in alloys). Therefore, in the 
following we will regard R c as a parameter ranging within 
given bounds for each element of interest. The other inde- 
pendent parameter of the model is the electron spacing r s , 
defined as 47rr|?/3 = N/V, where N/V is the conduction 
electron density. Such parameter enters Eq. J2J through 
the Fermi momentum kp in the Thomas-Fermi parame- 
ter /c 2 , and can be used as a measure of pressure, with r s 
decreasing as pressure increases. Moreover, for a given lat- 
tice structure, all inter-site distances Rij in Eq. JJJ scale 
with r s . 

As a result of the singular behaviour of h(q) dXq = 2kp, 
which in turn arises from the logarithmic discontinuity in 
the derivative of the Lindhard function x(q)i the Fourier 
transform in Eq. J2J endows the pair potential <£(i?) with 
an oscillating behaviour, with a characteristic length ~ 
ir/kp (Friedel oscillations, see Fig.^ inset). The analytical 
properties of the pair potential <P(R) have been analyzed 
further in Ref. fTS]. The actual value of the total energy 
U in Eq. (TjQ) then depends on whether the distances of 
nearest and farther neighbours in the lattice lay close to 
maxima or minima in the plot of <P(R) (Fig.rj). For a fixed 
value of the pseudopotential parameter R c , such distances 
can be actually shifted to lower values by decreasing r s , 
i.e. by means of an applied pressure. On the basis of such 
considerations, the stability of the crystal structures of 
several elements under pressure has been discussed in the 
past [231 



3 The case of lithium and other light alkali 
metals under high pressure 

The phase diagram of the alkali at low pressure has been 
extensively addressed both experimentally and theoreti- 
cally, with lithium crystallizing in the 9R phase at zero 
temperature and pressure 22 ,25 26 27 . At high pressure 
(~ 39 GPa), a new structural phase (Pearson symbol 
c/16) has been recently detected [7]. Hanfland et al. 9 
have recently observed bec — > fee — > c/I6 transitions in 
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Fig. 1. Scaled pair potential R<P(R) versus interatomic dis- 
tance R. Main plot: R c — 1.76 a.u., r a — 1.5 — 3.326 a.u., 
with thicker line corresponding to r s = 3.326 a.u. Symbols 
correspond to nearest and farther neighbours locations in the 
bcc (open circles), fee (full circles), and hep (triangles) lat- 
tices, respectively. Inset: Same as main plot for R c = 1.76 a.u., 
r a = 1.5 a.u., but now with distances scaled by k-p/n. Notice 
that minima in R<E>(R) occur close to integer values of kpR/ir 
(Friedel oscillations). 



Na. Thus we can regard the c/16 phase as our starting 
point for the following discussion on alkali metals. Such 
a high-pressure phase is characterized by a bcc primi- 
tive cell, with an 8-atom basis an d can be thought 
of as a distorted bcc phase, with a distortion parameter 
X — 0.045 — 0.060 0. The c/16 phase formally reduces to 
the bcc structure ('supercell' with eight usual cubic cells) 
for x = [TT]. In the undistorted bcc phase, the lattice 
is composed of two interpenetrating cubic sublattices, A 
and B, say. Such classification applies to the c/16 struc- 
ture as well, after distortion from the parent bcc lattice, 
but now with a basis of four atoms for each sublatticc. 
The c/16 phase by itself is not dimerized, and it has been 
recently predicted to be even superconducting 28,29 . We 
would like to test its stability with respect to a dimerized 
phase obtained by rigidly shifting the two sublattices each 
other of a tiny amount y along the (111) direction. Here, 
y plays the role of the order parameter discussed above. 
For a fixed value of the distortion parameter x, one can 
then think of expanding the total energy per atom Eq. 0} 
in powers of our 'dimerization' parameter y, 
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Fig. 2. Shaded regions in the plane of parameters (R c , r s ) are 
characterized by the condition U yy [x,0] < 0, which signals an 
instability towards dimerization. The regions corresponding to 
instability towards dimerization in the undistorted bcc phase 
(x = 0) and in the c/16 phase (x — 0.05) are bounded by con- 
tinuous and broken lines, respectively. Data for (R c ,r 3 ) pairs 
for the light alkali metals at ambient pressure have been taken 
from Refs. |24I23| . Different methods of fitting pseudopoten- 
tials to the experimental data result in a rather wide range for 
R c , which is here displayed as an error bar. 
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where n labels lattice points in the primitive bcc lattice, 
with lattice parameter a, /i and v label the basis vectors in 
sublattice A and B, respectively and iJjLj, denotes 
the mutual distances between lattice points for y = 0. 
Eq. (0 has been evaluated on a finite lattice, large enough 
to reach full convergence. Figure [5] displays our numerical 
results for U yy [x, 0] as a function of parameters (R c ,r s ), 
for x = (undistorted bcc phase) and x = 0.05 (represen- 
tative value of the high-pressure c/16 phase according to 
Refs. ITCH)- 

In the plane of parameters (R c ,r s ) we find regions 
U[x,y] = Uq + Uy[x,0}y + ^U yy [x,0}y 2 + Oiy 6 ), (5) f instability towards dimerization where U yy [x, 0] <0 

(Fig. shaded areas) . It is remarkable that at ambient 
pressure all the alkali metals are predicted to be stable 
in the symmetric phase. As the core radius R c is related 
to the band gap, one would in general expect R c to be a 
(monotonic) function of the density parameter r s . More 
generally, a more quantitative analysis would require a 
non-local pseudopotential, with a density-dependent range 
in reciprocal space. However, within linear response the- 
ory, the long wavelength behaviour of the pseudopotential 
is what actually matters. Indeed, the first zero of Eq. 0] 
occurs at q — ir/2R c , which is very close to the smallest 



where U y — dU /dy etc. Due to crystal symmetry, it can 
be proved analytically that U y [0, 0] = exactly, whereas 
we numerically checked that \U y [x, 0]/U yy [x, 0]| <C 1, for 
i<1 [s23- Therefore, an indication of instability towards 
'dimerization' is provided by the condition U yy [x, 0] < 0. 
From Eqs. and (J2J, U yy [x,0] can be expressed as 
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reciprocal lattice vector of most metals. Therefore, at the 
level of approximation implied by the present calculation, 
we can safely neglect the density dependence of R c , with- 
out qualitatively alter our main conclusions. In this way, 
the effect of applied pressure on the phase point corre- 
sponding to each alkali metal in Fig.[3]can be tracked as a 
vertical shift at constant R c , with increasing density. How- 
ever, the error bars for R c under normal conditions grossly 
provide order of magnitude boundaries for the generally 
expected variation of R c with density. 

A sign change in U yy for a given alkali metal then 
corresponds to a critical density. The uncertainty on R c 
lends a wide window for the critical densities of Li and 
Na: the largest values (still compatible with the core ra- 
dius data of Refs. |24U2H| ) are r s ~ 2 — 2.2 a.u., which are 
not too far from the values observed by Hanfland et al. 
|7j for the onset of a dimerized phase for Li (r s « 2.1 a.u. 
and P = 165 GPa) and predicted by Christensen and 
Novikov for the onset of a dimerized phase for Na 
(r s w 2.3 a.u. and P — 220 GPa). However, we caution 
that such a critical density occurs for r s ~ R c , i.e. close 
to the limits of applicability of the empty-core pseudopo- 
tential approximation, thus suggesting that more refined 
pseudopotentials should be employed in more realistic cal- 
culations. The phase diagram of Fig.[2]also suggests that, 
with a further increase of density, both Li and Na should 
drop below the instability area, and a second transition 
would eventually restore the symmetry of the lattice. This 
oscillating behaviour, from a symmetric phase to a less 
symmetric one and back to a restored symmetric phase, 
seems to be a general effect due to the presence of Friedel 
oscillations in the interatomic distance dependence of the 
pair potential. Such oscillations in turn are a consequence 
of the existence of a Fermi surface at k = fcp. Thus the 
tendency towards the formation of atomic pairs would be 
a signature of the Fermi liquid behaviour of the alkali met- 
als. 

Our method relies on linear response theory, which 
contains the Thomas-Fermi approximation as a long wave- 
length limit |31j . Both such methods are known to be re- 
liable if the density gradients are not too strong. In large 
atoms, the Thomas-Fermi method yields results which 
only differ by 8% from Hartree-Fock calculations, in space 
regions where the density gradient Vp/p « 2.5/eto (do is 
the Bohr radius) [H2j- Here, the typical density gradients 
are Vp/p < 0.5/a , as reported by Ref. ^ for the high 
pressure phase of Li at 165 GPa. Thus, we expect linear 
response to be reliable in the high pressure range, where 
eventual second order corrections are very small (for these 
gradients even the simple Thomas-Fermi approximation 
would deviate less than 1% from Hartree-Fock calcula- 
tions). 

4 Conclusions 

In summary, we have shown that a distorted bcc lattice 
(a/16 phase), surrounded by a Fermi liquid, undergoes 
several structural transitions under high pressure, oscil- 
lating between a symmetric phase and a broken-symmetry 



dimerized phase. These phase oscillations seem to be a di- 
rect consequence of the non monotonic behaviour of the 
pair potential which is characterized by a Friedel wave- 
length 7r/fcp in the presence of a sharp Fermi sphere. Such 
findings could be relevant for the understanding of the 
tendency towards the formation of atomic pairs, which has 
been recently reported for the light alkali metals. We fur- 
thermore predict that a restoring of symmetry should take 
place at some stage under higher pressure, which would 
be signalled by a reentrant metallic character. 

We thank N. W. Ashcroft, P. Ballone, N. H. March, G. Pic- 
citto, P. S. Riseborough, K. Syassen for useful discussions and 
correspondence. 
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